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ABSTRACT 

The lifetime of isolated protoplanetary disks is thought to be set by the com- 
bination of viscous accretion and photoevaporation driven by stellar high-energy 
photons. Observational evidence for magnetospheric accretion in young sun-like 
stars is robust. Here we report the first observational evidence for disk photoe- 
vaporation driven by the central star. We acquired high-resolution (R~30,000) 
spectra of the [Ne ll] 12.81//m hne from 7 circumstellar disks using VISIR on 
Melipal/VLT. We show that the 3 transition disks in the sample all have [Ne ll] 
line profiles consistent with those predicted by a photoevaporativc flow driven by 
stellar extreme UV photons. The ~6km/s blue-shift of the line from the almost 
face-on disk of TW Hya is clearly inconsistent with emission from a static disk 
atmosphere and convincingly points to the presence of a photoevaporative wind. 
We do not detect any [Ne ll] line close to the stellar velocity from the sample 
of classical optically thick (non-transition) disks. We conclude that most of the 
spectrally unresolved [Ne ll] emission in these less evolved systems arises from 
jets/outfiows rather than from the disk. The pattern of the [Ne ll] detections and 
non-detections suggests that extreme UV-driven photoevaporation starts only at 
a later stage in the disk evolution. 

Subject headings: accretion, accretion disks - infrared: stars - planetary systems: 
protoplanetary disks - stars: individual (TW Hya, CS Cha, T Cha, VW Cha, Sz 73, 
Sz 102, HD 34700) 
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Introduction 



It is well established that most ~1 Myr-old stars are surrounded by relatively massive 



optically thick dust disks. By an age o. 
retain an optically thick dust disk (e.g. 



■'lO Myr only a few 



Hernandez et al. 



percent of sun-like stars still 



20071 ) and many intermed iate-mass 



stars might already harbor second-generation dust disks (e.g. 



Currie et al. 



2008|). This 



fast clearing of primordial dust is in agreement with th e short formation timescales of 



chondrules, asteroids, and planets in the Solar System (Pascucci fc Tachibana 



200% . 



Although much less is known about the evolution of the gas component, there are at 
least three observables pointing to a similarly fast dispersal timescale: a) an order of 
magnitude lower accretion rate f or the few stars still accreting at an age o f ~10Myr in 



comparison to ~1 Myr-old stars (IMuzeroUe et al. 



2000 



Lawson et al. 



uppe r limits of less than .1 Mjnp in disks around non- accreting stars (IHollenbach et al. 



2005 



Pascucci et al. 



2004): b) gas mass 



20061 ) ; and c) upper limits on the H^-to-dust ratio of less than 10 in 



two ~ 12 Myr-old edge-on disks (iLecavelier des Etangs et al 



2001 



Roberge et al 



20051). 



Which are the physical mechanisms clearing out primordial disks? 

Models of protoplanetary disk evolution suggest that most of the disk mass is cleared 
out by viscous evolution ( accretion of gas on t o the central star) and photoevaporation driv en 



by th e central star (e.g.. 



Clarke et al. 



2001 



Alexander et al. 



2006b 



Gorti fc HoUenbach 



20091 ). The strongest evidence that young stars are accreting nebular gas comes from the 
veiling of optical and UV photospheric absorption lines w hich is produced by the continuum 



emission of an accretion shock at the stellar surface (e.g. 



Calvet et al 



2000). In addition. 



the large (few hundred km/s) line widths and asymmetries of permitted emission lines are 



found to be wel 



( iMuzeroUe et al. 



repro duced by models with infalling gas via magnetospheric accretion 



19981 ). In contrast to the numerous diagnostics of accretion, unambiguous 



diagnostics of centrally-driven disk photoevaporation are lacking. 
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Here we show that the profiles of [Ne ll] forbidden lines at 12.81 fim from disks with 
dust inner holes (hereafter, transition disks) are entirely consistent with those predicted 
from photoevaporative disk winds (Sect. Hj). These observations provide the first strong 
evidence in favor of centrally-driven disk photoevaporation. We also demonstrate that in 
disks with radially continuous optically thick dust (classical, non-transition disks) most 
of the [Ne ll] emission detected with the Spitzer Space Telescope does not arise in a disk 
but likely in a jet/outflow (Sect. [5]). Finally, we discuss how observations, such as those 
presented in this paper, can be used to identify the onset of disk photoevaporation (Sect. [6]) 



2. Motivation, Observations, and Data Reduction 



The IRS spectrograph on board the Spitzer Space Telescope has recently detected 
the [Ne ll] 12.81 /xm line toward many young protoplanetary disks as well as toward the 



possibly more evolved transition disks ( 



2007 



Espaillat et al. 



2007 



Geers et al. 



Ratzka et al. 



2007 



2006 



Pascucci et al 



Flaccomio et al. 



2OO9I : 



2007 



Lahuis et al 



Guedel et al 



20091). 



Lines appear spatially and spectrally unresolved indicating that the [Ne ll] emission 
is confined to within ~1,000AU from the central star and lines are narrower than 
~500km/s. Models of disks irradiated by stellar X-rays or extreme UV (EUV, hu > 
3.6 eV) photons can sufficiently ionize the disk surface to reproduce the observed line fluxes 



(jGlassgold et al 



2007 



HoUenbach fc Gorti 



Meiierink et al. 



2008 



Gorti fc HoUenbac 



2 



2008 



Ercolano et al 



2008 



20091 ) suggesting that [Ne ll] emission is a robust tracer of disk gas. 
However, th e flrst two spectrally resolved [Ne ll] lines brought some surprises. In the case 



of TW Hya, 



Herczeg et al. 



20071 ) found that the [Ne ll] line is centered at the stellar radial 



velocity, consistent with a disk origin, b ut about a factor o: 



from X-ray and EUV irradiated disks. 



van Boekel et al. 



two broader than that predicted 



(I2OO9I ) spatially and spectrally 



resolved the strong [Ne ll] line from the T Tau complex demonstrating that in this system 
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most of the [Ne ll] emission originates in jets/outflows with neon atoms partly ionized by 
stellar X-rays. 

Our aim here is to enlarge the sample of spectrally resolved [Ne ll] lines and understand 
in which systems [Ne ll] emission originates in a disk. In the following subsections 
we describe our observational campaign and the data reduction (Sects. 12.11 12.21 12.31) . 
Observations were carried out with the high- resolution (R~30,000) spectrograph VISIR 
mounted on the VLT telescope Melipal (Sect. 12.21 for more details). In addition to the 
already published Spitzer/IRS spectra, we reduced and present here new archival IRS 
spectra that aid the interpretation of the VISIR observations (Sect. 12.31) . 



2.1. Target Selection 



Because of the lower- sensitivity of VISIR with respect to Spitzer/IRS, we restricted 
our sample to disks with bright/unresolved [Ne ll] lines (Spitzer fluxes greater than 
10^^^ ergcm"^ s^^). In addition, we selected disks brighter than ~50mJy in the continuum 
because target acquisition for fainter objects becomes challenging with VISIR. The 6 targets 
we selected for this campaign are TW Hya, CS Cha, VW Cha, T Cha, Sz 73, and Sz 102 
(see Table [1] and IHfor their main properties). At the beginning of the second night we 
also observed the double-lined spectroscopic bi nary HP 34700A, w hich is the brightest 



infrared T Tauri couple in a quadruple system (ISterzik et al. 



20051 ). Although HD 34700A 



did not have a published [Ne ll] line detecti on, its near- and m id-infrared spectra present 



unusually strong PAH e mission bands (e. 



stellar far-UV flux (e.g. 



Geers et al 



Smith et al. 



20041 ). possibly hinting to a high 



20061 ). which may be indicative of a high EUV flux. 



The systems we selected can be grouped into three categories based on their broad-band 
spectral energy distribution (SED). VW Cha, Sz 73, and Sz 102 have significant excess 
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emission relative to the photospheric flux f r om near- through to far-in frared wavelengths 



( IGauvin fc Strom 



1992 



Hughes et al. 



199 



Kessler-Silacci et al 



20061). Such broad SEDs 



can be well reproduced by radially continuous optically thick dust disks, extending from 
a few stellar radii out to hundreds of AU. TW Hya, CS Cha, and T Cha are classified 
as transition disks in the literature. Their SEDs present a strongly reduced (or lack of) 
near-infrared excess emission but large mid- and far-infrared emission. Detailed modeling 
of their SEDs points to relatively large inner dust cavities almost devoid of sub-micron- and 



micron-sized dust gr ains: ^1-4 AU for TW Hya flCalvet et al. 



-43 AU for CS Cha (lEspaillat et al. 



2002 : 



jlatzka et al 



20071), and ~15AU for T Cha (IBrown et al 



2007) 



2003). 



There is no evidence of gas in ner holes in these syst ems: TW Hya is accreting disk gas at 
a rate of ~ 5 x 10~^° Mr^/vi (|Muzerolle et al.ll2000l ). the spectroscopic binary CS Cha at 



< 10-s Mq/yv flEspaillat et al. 



2007 



while the small (< lOA) but v ariable Ha equivale nt 



width from T Cha suggests low-level and possibly episodic accretion (lAlcala et al 



19931 ). 



Finally, the SED of HD 34700A has little excess emission at a 



tenuous dust disk, possibly a debris disk (jSylvester fc Skinner 



1 wav elengths produced by a 



19961 ). 



2.2. Observations 



We performed long-slit high-resol ution spectroscopy with the spectrograph VISIR 



mounted on the VLT telescope Melipal fiLagage et al. 



20041 ). The observations were executed 



on 20, 21, and 22 February 2008. VISIR has a 256 x 256 BIB detector array with an image 
scale of O'.'127/pixel. We acquired all targets with the PAH2-NEII filter. We used a slit 



width of 0'.'4 centered at 12.81/im to detect and spectrally resolve the [Ne ll] (^P-- 



3/2 P1/2) 



4he estimate from 



Espaillat et al. 



(I2OO7I ) should be considered only an upper limit be- 



cause the UV excess emission was computed assuming that CS Cha is a single star 
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fine-structure line at 12.81355/xm (lYamada et al.lll985l ). With this configuration we covered 



the spectral region between ~ 12.79 and 



measured from first light data (iKaufl 



-2.83 /im with a resolution of R ^ 30, 000 as 
20061 ) and from the FWHM of 5 narrow O3 sky lines 
from our spectra. This resolution corresponds to 3 VISIR pixels or ~10km/s in velocity 
scale. 

We applied the standard chopping/nodding technique to suppress the mid- infrared 
background. For each target the slit was positioned in the default North-South orientation 
and chopping/nodding was done along the slit with a throw of 8". Total on-source exposure 
times per science target were 1 h or longer. To correct the spectra for telluric absorption 
and to obtain an absolute fiux calibration, we observed standard stars immediately before 
or after the science target and at a similar airmass. A summary of the observations is 
presented in Tabled 



2.3. Data Reduction 

VISIR high-resolution spectra. For each nod position VISIR raw data are 
organized in cubes with the third axis having 2n + 1 planes, where n is the number of 
chopping cycles. For each chopping cycle (i) two so-called Half-Cycle exposures are made: 
the Ai image from the on-source position of the chopper, and the Bi image from the 
off-source position of the chopper. The data cube stores the Ai images in the odd planes 
and the average of the current and all previous Ai — Bi images in the even planes. The last 
plane of the cube contains the average of all Ai — Bi images and it is thus identical to the 
2 X n plane. 



We process the raw data cubes with the VISIR pipeline version 3.2.1 (Lundin et al 



20081 ) . First, 2D frames in each nod position are created by averaging all the Half-Cycle 
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difference images in the data cube and low-level vertical stripes are removed. This step 
results in a set of images with a positive and a negative beam because the chop throw 
we chose is smaller than half of the VISIR field-of-view. Then, nodded images (images 
corrected for the telescope nodding) are produced by averaging the images in the two nod 
positions (A and B) and dividing them by 2 times the detector integration time. At this 
point each nodded image contains a double positive beam in the center and a negative 
beam on each side separated by the nodding(=chopping) throw. Pixels detected as bad in 
the Half-Cycle frames are now cleaned by interpolation with neighboring pixels. In parallel, 
a reference frame of the infrared background, not corrected for chopping or nodding, is 
also created from the first Half-Cycle of the first nod position. The nodded images as well 
as the infrared background image are corrected for the optical distortion, which is known 
analytically, and are shifted and added to form the final combined image using the offsets 
stored in the FITS header (the nodding cycle sequence used for VISR is AnBnBnAn). The 
last step before extraction is to fold the two negative beams into the central positive beam. 

The spectrum is then extracted following the optimal extraction method by [Horns 



( 1l986l ). In brief, the image is collapsed in the spectral dispersion direction to create a source 
profile, which is normalized and then expanded in spectral direction to obtain a weight 
map. The science frame is multiplied by the weight map and the spectrum is obtained 
by summing in spatial direction. This extraction method works well for sources where, in 
addition to narrow emission line(s), the continuum emission is also detected. For sources 
like CS Cha, where only a narrow emission line is detected (but not the continuum, see 
Fig. [1]), the optimal extraction method fails. For this source and its calibrators we have 
performed the standard full aperture extraction (same weights within a defined spatial 
range of pixels). We tested that an aperture of 10 pixels gives the best S/N in the extracted 




spectra. 
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The wavelength cahbration is done using the collapsed background frame in spatial 
direction and cross-correlating it to a synthetic model spectrum of the atmosphere. The 
dispersion relation is well approximated by a first order polynomial and the offset in 
pixels that maximizes the cross-correlation can be determined to an accuracy of 0.01 pixel 
(~0.03km/s). We also compared the observed peak position of two photospheric absorption 
lines in the spectra of the K5 III star HD 136422 and the K4.5 III star HD 139127 to 
those from the MARCS model atmosphere of a K5 III star (HR 6705, the spectrum was 
kindly provided by L. Decin). This exercise demonstrates that the peak centroids can be 
determined to a lower accuracy, varying from a tenth up to a few km/s (see also Table 3). 

All steps described above have been applied to both the target and its standard star(s). 
In addition, the extracted spectrum of the standard star was flux calibrated using the model 
flux from the VISIR standard star catalog. From this calibration we created a spectral 
response function (Jy/ADU) that we applied at the wavelengths of the science spectrum. 
This step removes most of the "fringing" present in the VISIR spectra (the flux modulation 
with wavelengths has a peak to pe ak amplitude of about 15% and is found to be stable over 



periods of at least several months. 



van Boekel et al. 



20091). 



Spitzer/IRS spectra. As mentioned in Sect. 12.11 all targets except HD 34700 
have published low- (R~105) or high-resolution (R~600) IRS spectra with bright [Ne ll] 
emission lines (see Table [1]) . We have reduced the archival IRS low-resolution spectrum of 
HD 34700 and report no detection of the [Ne ll] line. In addition, we have reduced archival 
high-resolution IRS spectra f or TW Hya and CS Cha obtained at different epochs than the 



already published spectra by 



Ratzka et al. 



(120071 ) and 



Espaillat et al. 



(120071 ). These spectra 



were acquired as part of the Spitzer GO program 30300 (PI, J. Najita) and are presented 



here to aid the interpretation of the VISIR results (Sect. HI). The da ta reduction o: 



high- and low-resolution spectra followed the procedure outlined in 



Pascucci et al 



the IRS 



(12006 . 



20071 . |2008| ) and iBouwman et al.l (120081 ) . The [Ne ll] line fluxes from these new spectra are 
also reported in Tabled] (first entries in column 5). 



3. Results 

We detect and spectrally resolve [Ne ll] emission lines from 4 out of 7 targets, 
specifically TW Hya, CS Cha, T Cha, and Sz 73. The [Ne ll] emission detected with 
VISIR is found to be comparable to the angular resolution (~0.5") estimated from both 
the mid-infrared continuum of our targets and from the telluric standards observed before 
or after the targets. The profiles of the lines are consistent with a single Gaussian profile 
(see Fig. [2]). Table [3] provides the peak centroids (in the stellocentric frame), FWHMs, and 
fluxes of the [Ne ll] lines computed assuming a Gaussian profile and a first-order polynomial 
for the continuum. In the case of non-detections (see Fig. [3]), we fit a first-order polynomial 
within ±100km/s of the star velocity. Table [3] provides the 3cr upper limits to the [Ne ll] 
flux computed from the RMS in the baseline-subtracted spectrum and assuming a line 
width of lOkm/s, equal to the instrument resolution (Sect. 12. 2p . Upper limits for broader 
lines simply scale by (FWHM/10). Figs. [2] and [3] illustrate our best fits to the detected lines 
and the hypothetical 3a upper limits when no line is detected. 

The second-epoch Spitzer spectra of TW Hya and CS Cha presented here demonstrate 
that [Ne ll] line fluxes are not constant in time and that the infrared continuum emission 
is also variable. More details on the line and continuum variability of TW Hya will be 
presented in an upcoming paper (Najita et al in prep.). The important result for the 
interpretation of these VISIR spectra is that changes of ~30% in the line and/or in the 
continuum are possible even among the class of transition disks. 

We do not detect the [Ne ll] emission line in the Spitzer low-resolution spectrum of 
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HD 34700A. This star is surrounde d by a tenuous dust d isk and is thus more similar to the 



optically thin dust disks studied by 



Pascucci et al. 



(I2OO6I ) rather than to the protoplanetary 



disks presented in this paper. The [Ne ll] non-d etection suggests tha t very little or no 



20071 ). corroborating its 



primordial gas is left in the disk of HD 34700A (jPascucci et al. 
identification as a debris disk. 

In the following, we discuss in more detail the 6 protoplanetary disks with spectrally 
unresolved [Ne ll] lines in their Spitzer spectra. We separate transition objects from radially 
continuous dust disks because their VISIR spectra appear remarkably different. 



4. [Ne 11] emission from the disk of transition objects 



For the three transition disks TW Hya, CS Cha, and T Cha we detect and spectrally 
resolve [Ne ll] lines located near the stellar velocity. In the case of CS Cha we detect 
only the emission line but not the continuum^. For all sources [Ne ll] line fluxes from 
VISIR spectra are lower than those measured from Spitzer spectra of factors between 
~0.5-0.8. One possibility is that the [Ne ll] emission is extended beyond the 0.4" slit 



adopted in these observations. Spat ially extended [^ 



outflow sources, as demonstrated by 



van Boekel et al. 



e 11] em ission can certainly arise in 

and by our sample of optically 



thick d ust disks (see Sect. [5 l) . However, no jets have been reported toward TW Hya and 



T Cha (lAzevedo et al 



2007 



Alcala et al 



in Ha has been detected toward CS Cha ( Takami et al. 



1993) and on 



y a c ompact (8-10 mas) outflow 



2OO3I ) . In Sect. 14.11 we also show 



that jets/outflows are not likely to be the main source of [Ne ll] emission in transition 



^The flux calibrated spectrum shown in Fig. [2] is the mean of the flux calibrated spectra 
from the 3 datasets, errors at each wavelength are the standard deviations from the 3 spectra 
on a common wavelength scale 
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objects. Another explanation for the Spitzer-VLT flux difference is line and/or continuum 
variability. At least for CS Cha and TW Hya we know that line/continuum variability 
measured from 2-epochs of Spitzer spectra is ~30%, similar to the [Ne ll] flux differences 
from the VLT and Spitzer spectra. We should also keep in mind that at the distance of 
TW Hya and T Cha (Table Hj) our slit just covers out to ~10AU from the central stars. As 
discussed later even [Ne ll] emission from the disk surface can extend beyond several tens 
of AU. Further observations with wider slits and different orientation angles are necessary 
to constrain the spatial extension of the [Ne ll] emission. 

The most important result from the VISIR spectra is that the measured [Ne ll] lines 
of transition disks are relatively narrow (ranging from ~14km/s to ~40km/s) and peaked 
near the stellar velocity (see Table [3]). These line characteristics are consistent with a disk 



origin for the emiss ion 



Pontoppidan et al 



. We note that TW Hya, whose disk is almost seen face-on (4 ± 1°, 



20081 ). has the narrower line width. The inclinations of the disks aroun d 



CS c 



la and T Cha are not well constrained but likely closer to edge-on. 



Espaillat et al 



(120071 ) successfully model the SED of CS Cha with a disk viewed at an inclination of 



60°. T Cha has a high projectional velocity suggesting that it is viewed almost edge-on 



60°. T Cha 1 


las a 


(Alcala et al. 


1993) 



viewing angle: the narrowest [Ne ll] line is measured from the almost face-on disk of 
TW Hya, a broader line from the inclined disk of CS Cha, and the broadest line from the 
almost edge-on disk of T Cha. This trend also suggests that the [Ne ll] line profiles are 
dominated by Keplerian rotation rather than by turbulence. 

Further constraints on the region traced by t he [Ne ll] emission can come from modeling 



line profiles. In a recent study. 



Alexanded (120081 ) showed that the velocity structure of a 



^The line width from TW Hya reported by lHerczeg et al.l (120071 ) is about 1.6 times larger 
than our. Their spectrum has a poor S/N and was obtained at an airmass as large as 1.8. 
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photoevaporative disk wind results in a line profile that is different from that produced 
by a bound disk atmosphere. In particular, when a disk is viewed face-on the vertical 
component of the wind velocity leads to a profile that is about ~10km/s in width (broader 
than the thermal line width) and blueshifted by ~6-7km/s. When a disk is viewed closer 
to edge-on, line profiles become broader (~ 30 — 40km/s), are double peaked, dominated 
by Keplerian rotation, and centered at the stellar velocity. Fig. H] compares the observed 
[Ne ll] line pr ofiles with the line prof iles predict e d by the standard photoevaporative disk 



wind model of 



Font et al.l tooi ) and 



Alexanderl (120081 ). the one in which there is no hole in 



the gas disk. We have chosen this model because our transition objects are still accreting 
disk gas, suggesting that there is no hole in the gas disk even if the inner dust disk is 
largely depleted of small dust grains (Sect. [2]). The standard photoevaporative disk wind 
model assumes that the disk atmosphere is heated by the stellar EUV photons at 10^ K and 
requires just 3 input parameters to compute the [Ne ll] line profiles: the stellar mass, the 
disk inclination, and the stellar ionizing fiux ($). The profiles shown in Fig. [Hare for the 
star /disk parameters summarized in Table H] and assuming a stellar ionizing flux of 10^^ 
photons/s. There is an almost perfect match between the observed and predicted [Ne ll] 
profile from the disk of TW Hya. The almost face-on disk of TW Hya has a blueshifted 



[Ne ll] line fully consi stent with 



toward the observer ( lAlexander 



; he em ission originating in a photoevaporative wind moving 



20081 ). A static disk atmosphere would result in a line 
emission centered at the stellar velocity, clearly inconsistent with our observations. Thus, 
we consider this result as the first strong evidence for a disk being photoevaporated by its 
central star. 

The [Ne ll] line profiles of CS Cha and T Cha are broadly consistent with predictions 
from the photoevaporative disk wind model (Fig. Hj), but there are also a few discrepancies, 
especially for the spectroscopic binary CS Cha. For this binary the predicted [Ne ll] line is 
broader than observed. This could result from a disk viewed at less than 45° inclination. 
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but the emission should be shghtly more blue-shifted than observed. 



In the simplest assumption that neon atoms are ionized only by stellar EUV photons, 
we can use the observed [Ne ll] line luminosities to estimate the stellar ionizing flux $ 
reaching the disk surface. The reason is that [Ne ll] line fluxes scale approximately linearly 
with $ (for $ < 10^^ photo ns/s) while the shape of the line proflle is insensitive to the 



value of $ (lAlexander 



20081 ). The comparison between the observed and predicted line 
luminosities yields stellar ionizing fluxes of 2.5, 13, and 1.9x10^^ photons/s for TW Hya, 
CS Cha, and T Cha respectively. These values have an uncertainty of a factor of 2-3 
dominat ed by the uncerta inty in the fraction of neon atoms that exist as Ne"*" in the disk 



surface. 



Alexanderl (120081 ) keeps this fractior i constant to 1/3 w^ 



models suggest that it could be as high as 1 ( jHoUenbach fc Gorti 



lilemore detailed gas 



20091 ). The ionizing rates 



we compute from the [Ne ll] lines seem plausible for accreting T Tauri stars. For TW Hya 
far-UV observations suggest values of 7 x 10^^ or 5 x 10^^ depending on h ow much EUV 



2003) 



emission from the accretion shock is attenuated by the accretion streams (iHerczegl 
As a next step we can use the values of $, in combination with star and disk masses, to 
estimate the disk wind rates and the disk lifetimes if EUV-d riven photoevapor ation would 



Font et al. 



(12004h we flnd disk 



be the only disk dispersal mechanism. Using the results from 
wind rates ranging from ~ 2 to 5 x 10~^°M(7)/y r. At these rates only the relatively low-mass 



disk of T Cha (0.003 Mq, 



Lommen et al 



20081 ) could be photoevaporated in ~10 Myr. The 



disks of TW Hya and CS Cha would disperse in several tens of millions of years. However, a 
more realistic estimate for the disk lifetime should include on-going mass loss onto the star 
via viscous accretion. In addition, disk lifetimes can be further reduced once the viscous 
accretion infall rate falls below the photoevaporation rate. At this point a gap opens in the 
gas disk and the direct EUV flux from the star is expected to disperse the entire disk in 
only ~ 10^ yr. As we discuss in Sect. [6l at least two of our transition objects may be on the 
verge of opening an inner gap in the gas disk. 
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Finally, it is worth commenting on the extension of the [Ne ll] emission in the 
photoevaporative disk models. Calculations from Alexander (priv. comm.) show that more 
than 90% of the [Ne ll] emission arises inside ~ 2 x 9(M^/M0) AU. These radial distances 
are covered by our 0.4" slit for CS Cha, are comparable to half the slit width for TW Hya, 
and smaller by a factor of 2 for T Cha. This shows that we certainly missed some [Ne ll] 
flux (especially from the disks of TW Hya and T Cha) even in the case that ionization is 
solely fro m EUV photons. Stella r X-rays could contribute to the gas ionization at larger 



disk radii (iGlassgold et al 



photoevaporative flow (lErcolano et al. 



2007) and t 



le X- ray heated gas might also participate in the 



20091 ). If X-rays substantially contribute to the 



[Ne ll] fluxes then the stellar ionizing fluxes we calculated above should be viewed only 
as upper limits. It would be extremely interesting in the future to compare the observed 
[Ne ll] profiles to those predicted by photoevaporating X-ray-heated gas disks. 



4.1. Can jets explain the [Ne ll] emission from transition disks? 



van Boekel et al. 



(120091 ) demonstrate 



The observations of the T Tau complex from 
that [Ne ll] emission can also originate from shocks generated by protostellar outflows/jets. 
Because of the high ionization potential of neon atoms (21.6 eV), substantial ionization can 
be only produced by fast shoc ks with speeds >70km/s f or typical pre-shock densities of 



10^ cm ^ (also called J shocks, iHoUenbach fc McKee 



198 



^or this type of st r ong r adiative 



van Boekel et al. 



(120091 ) argue 



shocks the shocked gas moves at almost the shock velocity, 
that these high velocities may be reached in the outflow from the T Tau S source which has 
a sight velocity of only ~40km/s, but it is likely to be in the plane of the sky. We now show 
that it is very unlikely that jets/outflows could explain the [Ne ll] emission of transition 
disks. First, let us consider the case of TW Hya. Its disk is seen almost face-on so we 



would expect any jet/outflow to be almost in the direction of the observer and the sight 
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velocity to be close to the actual shock velocity. The sight velocity of the [Ne ll] emission 
from TW Hya is just -6km/s and the line is relatively narrow (deconvolved FWHM of 
~10km/s). If this velocity is representative for the shock velocity as predicted from J shock 
models, it would be unabl e to appreciably ionize neon atoms. Corroborating this assertion, 



HoUenbach &: Gortil (l2009l ) show that the [Ne ll] line luminosity produced by the postshock 
region of a radiative shock depends linearly on: M^v"^, where My^ is the protostellar mass 
wind loss rate and Vg is the shock velocity. The wind mass loss rate s cales with the stellar 



mass accretion rate as ~ 0.01 — 0.1 x Mace (IHartigan et al. 



1995 : 



White &: Hillenbrand 



20041 ). Given the low stellar accretion rates of transition disks, the [Ne ll] fluxes reported in 
Table [3] could be reached only with shock velocities ranging from several hundreds of km/s 
to thousands of km/s, clearly inconsistent with the small velocity shifts and relatively small 
FWHM of the observed [Ne ll] lines. Finally, line widths from shocked gas are expected 
to be quite broad (s everal tens of km/s) regardless of the observer viewing angle (e.g.. 



Hartigan et al 



19871 ). thus shocked gas cannot reproduce the observed trend between the 



FWHM of [Ne ll] lines and the disk inclination. Although shocks generated by jets/outflows 
are unlikely to explain the [Ne ll] emission from transition dis ks, they may domin ate the 



[Ne ll] emission from less-evolved optically thick dust disks (see 
Section). 



Guedel et al. 



2009 



and next 



5. Jets/outflows dominate the [Ne ll] emission of classical optically thick disks 

The three optically thick dust disks observed here have spectrally unresolved [Ne ll] 
fluxes comparable to those of transition disks. Nevertheless, we do not detect any [Ne ll] 
line close to the stellar velocity in their VISIR spectra. Only in the case of Sz 73 we detect 
a broad (~60km/s) emission line at ~-100km/s with respect to the stellar velocity and 
a line flux that is just ~30% lower than the Spitzer [Ne ll] flux. The large FWHM in 
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combination with the large blueshift point to a jet/outflow origin for the [Ne ll] emission. 
The case of Sz 73 well illustrates that when jets/outflows are present the [Ne ll] emission 
they produce can domin ate over the disk emissio n. This is in line with the observations of 
the T Tau complex from 



van Boekel et al. 



(120091). 



The [Ne ll] non-detections of Sz 102 and VW Cha are also consistent with this scenario. 
If the [Ne ll] emission measured with Spitzer would originate in a disk, only relatively broad 
lines (FWHM >50km/s), compatible only with close to edge-on disks, would be undetected 
in our VISIR spectra. Interestingly, VW Cha and Sz 102 are k nown to have out f lows at 



posit i on angles (from N to E ) of 90 and 99 degrees respectively (iBally et al. 



19861: 



Graham fc Heyer 



2006 



Krautter 



19881 ). Since we did not change the default slit orientation of VISIR 
(PA=0°) we acquired spectra in a direction perpendicular to the known outflows. Thus, 
our non-detections can be reconciled with the Spitzer observations if most of the [Ne ll] 
emission arises in jets/outfiows, which can be verified with future high-resolution spectra 
covering the region of the known outflows. 



6. On the onset of centrally-driven photoevaporation 



Viscous accretion is thought to be the dominant disk dispersal mechanism but alone 



cannot explain the relatively short (< lOMyr 



observations (e.g. review by 



DuUemond et al. 



disk dispersal timescales inferred from 



20071 ). High-energy stellar photons can 



heat and ionize the disk surface and induce photoevaporation (i.e. disk mass loss) from 
outer disk regions. Exactly when photoevaporation becomes an important disk dispersal 
mechanism and how much disk mass it can remove is still a matter of debate. FUV 
and X-ray photons can penetrate much larger column densities of gas than stellar EUV 
photons. As a consequence, they are expected to drive efficient photoevaporation from the 
early times of disk evolution at rates that may exceed the EUV photoevaporation rate by 
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one or even two orders of magnitudes (iGorti fc HoUenbachl l2009l : lErcolano et al.l |2009| ) . 
Stellar EUV photons start to penetrate the disk wind, and thus reach the disk surface, 
only when accretion rates fall below ^ 10~^Mf7^/yr, correspon ding to hydrogen column 



density screens of N(H)<10^^cm ^ (IHollenbach fc Gorti 



200911. Although EUV photons 



20091 ). they can help to 



likely produce little mass loss at all stages (iGorti &: HoUenbachl 
quickly (~ 1 0^ yr) disperse remaining disk gas once stellar accretion rates fall below a few 



lO-^°M0/yr flAlexander et al. 



2006ah . 



Recently, 



Currie et al 



(120091 ) have argued that the large number of evolved dust disks 
in the 5 Myr-old cluster NGC 2362 is inconsistent with predictions from the standard UV 
photoevaporation model, in which EUV-dri ven photoevaporation follows viscous accretion 



( IClarke et al. 



2001 



Alexander et al 



2006bl ). Indeed, these models predict a paucity of 
evolved disks because once the accretion inflow rates fall below the photoevaporation rate 
the disk dispersal is fast (~ 10^ yr). However, model predictions are directly relevant only 



to the evolutioi i of the gas disk 



dust evolution, Alexander et al. 



compon ent. In order to translate the gas evolution into 



(j2006bl ) made the simplistic assumption that the dust 



is coupled to the gas throughout the entire disk evolution. This is cer tainly unrea^ 



because grain growth, observed in almost all protoplanetary disks (e.g. 



Natta et al. 



istic 



2007h. 



results in dust settling, i.e. decoupling of the dust from the gas disk ( iDuUemond fc Dominik 
2OO4I ). A more direct way to test photoevaporation models is to constrain the evolution of 



the gas disk through observations of gas lines tracing the accreting as well as the possibly 
photoevaporating material. In what follows we will take this approach to speculate on the 
onset of photoevaporation. 

We start by discussing our sample of probably less evolved classical optically thick 



disks with radially continuous dust distribution. The Ha equivalent widths 



three systems are greater than 50 A ( Hughes et al 



1994 



Guenther et al 



rom these 



20071 ). indicating 
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that the central stars are accreting nebular gas. This is corroborated by the detection of 
outflo ws and large near-inf rared excess emission, both characteristic to actively accreting 
disks (IHartigan et al.l 119951 ). In these systems we do not find evidence for on-going disk 
photoevaporation. The [Ne ll] line flux upper limits at the stellar velocity imply that 
not only EUV photons but also stellar X-rays are efficiently screened by the protostellar 
wind and accretion colu mns, requiring large gas column densities {Nh > lO^^cm"^, 



e.g. 



Guedel et al 



20071 ). If column densities exceed ~ 10 cm , FUV photons are 



also significantly attenuated (IHoUenbach fc Gorti 



20091 ). For lower column densities. 



photoevaporation driven by FUV photons could be occurring and we would not detect 
it because FUV photons have energies lower than the ionization potential of neon atoms 
(21.6 eV). These systems may be at the same evolutionary stage as the T Tau complex. 



Recently, 



Najita et al. 



(120091 ) reported spectrally resolved [Ne ll] lines from two disks 
in the Taurus- Auriga star-forming region: AA Tau, a classical optically thick disk whose 



star presents periodic accretion rate bursts (iBouvier et al. 



20071 ): and GM Aur, a transition 



disk whose c entral star is accreting at a level comparable to that of classical T Tauri stars 



;iO-SM0/yr, 



GuUbring et al. 



19981 ). The [Ne ll] lines are found to be centered at the stellar 
velocity suggesting that the emission originates from a static disk atmosphere very likely 
heated by stellar X-rays. However, the measured fluxes are factors of 2 and 3 lower than the 
spectrally unresolved fluxes from Spitzer spectra. One possible explanation is that most of 
the undetected [Ne ll] emission originates far from the central stars, in ambient gas shocked 



by a jet/outflow. This seems very plausible for AA Tau, whose optic al spectrum 
forbidden emission lines typical of high- velocity jets (e.g. [Nil] 6530 A, 



nas 



Hirth et al. 



19971) 



and a jet has als o been imaged with HST/STIS and Goddard Fabry Pe rot spectrograph 



Najita et al. 



(120091 ) does not seem 



( ICox et al.ll2005l ). The [Ne ll] emission reported by 
to trace photoevaporating gas in these two systems, possibly because X-ray heated gas 
producing [Ne ll] emission is slightly cooler (~5,000K) and thus more bound than the 
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EUV-heated gas. 

Finally, the transition disks presented here have [Ne ll] line profiles consistent with 
emission from a photoevaporative disk wind indicating on-going photoevaporation at this 
evolutionary stage. In the case of TW Hya, the stellar accretion rate is close to the limit 
where EUV-driven photoevaporation could start opening a gap in the gas disk. The small 
but varying Ha equivalent width of T Cha suggests that this transition object may also 



be on the verge of clearing a gas hole. Interestingly, TW 



to be relatively old, older than 5 Myr fiWebb et al. 



1999 



Hya and T Cha are a 



van den Ancker et al. 



so th ought 
1998, see 



also Tabled]). CS Cha is the youngest among our transition disks. Its upper limit to the 
stellar accretion rate is close to what is needed for EUV photons to start penetrating the 
disk wind, placing it at an earlier evolutionary stage than TW Hya and T Cha (though the 
evolution of the inner disk is certainly affected by presence of the stellar companion). 

In overall, these results support the evolutionary picture in which EUV-driven 
photoevaporation occurs at later times in the disk evolution and the mass loss rates can 
be only modest. We should also note that if stellar X-rays substantially contribute to the 
[Ne ll] emission of transition disks then the stellar ionizing flux $ calculated in Sect. H] 
would be only upper limits, further reducing the efficiency of EUV-driven photoevaporation. 
Future observations should aim at expanding the sample of classical and transition disks 
with spectrally resolved [Ne ll] lines and accretion rate measurements. In addition, it would 
be extremely valuable to compute [Ne ll] line profiles from X-ray heated photoevaporating 
gas that could be directly compared to observations. 



7. Summciry 



We used the high-resolution spectrograph VISIR on the VLT to spectrally resolve 
[Ne ll] emission hues previously detected with the Spitzer Space Telescope in a sample of 
young protoplanetary disks. Our main results can be summarized as follows: 

• We detect and spectrally resolve the [Ne ll] hues from the 3 transition disks in our 
sample. The hne centroids, widths, and intensities indicate that the [Ne ll] emission 
arises in a disk. Models of disk wind photoevaporation well reproduce the observed 
line profiles strongly suggesting on-going photoevaporation driven by stellar EUV 
photons at this evolutionary stage. 

• Prom the sample of optically thick radially continuous dust disks we only detect a 
strong and broad [Ne ll] emission blueshifted by -lOOkm/s toward Sz 73. This points 
to a jet/outflow origin for the [Ne II] emission. The [Ne II] non-detection in the other 
two systems with known jets also suggest that most of the spectrally unresolved 
[Ne ll] emission detected with Spitzer does not arise from the disk but rather from 
surrounding gas shocked by jets/outflows. 

These results provide the first observational evidence for centrally-driven disk photoe- 
vaporation. Because transition disks have already undertaken significant evolution in 
comparison to the radially continuous dust disks, our results demonstrate that EUV-driven 
photoevaporation can occur only at a late stage in the disk evolution. 

We thank R. D. Alexander for proving the [Ne ll] line profiles from photoevaporative 
disk winds and for extremely useful discussions. We would also like to thank F. Lahuis 

for making available the high-resolution Spitzer/IRS spectra of Sz 73, Sz 102, T Cha, and 
VW Cha and L. Decin for the MARCS model atmosphere of HR 6705. We are also grateful 
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to D. HoUenbach, U. Gorti, B. Ercolano, J. Najita, J. Carr, and D. Apai for valuable 
discussions, and an anonymous referee for helpful suggestions. IP is pleased to acknowledge 
support through NASA contract 90035375. 
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Fig. 1. — Combined image of CS Cha from the first dataset (February 22). The [Ne ll] 
emission is clearly detected both in the positive and in the negative beams but there is no 
detection of the continuum emission. The combined images from the other two datasets are 
similar and therefore not shown here. 
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Fig. 2. — Sources for which we detect [Ne ll] emission hnes in their VISIR spectra. The 
X-axis gives the velocity in the stellocentric frame (see Table [T] for the stellar heliocentric 
radial velocity). On top of the continuum emission we overplot the best Gaussian fits to the 
data (dashed lines). In the case of CS Cha the continuum emission is not detected (see also 
Fig.T). 
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Fig. 3. — Sources for which we do not detect any [Ne ll] emission hne in their VISIR spectra. 
The X-axis gives the velocity in the stellocentric frame. On top of the continuum emission 
we overplot the hypothetical 3a upper limits (dashed hues) reported in Table [3] for line 
widths equal to lOkm/s. Note that spectrally unresolved [Ne ll] lines have been detected 
with Spitzer toward VW Cha and Sz 102 (Table [1]). Our non-detections suggest that most 
of the [Ne ll] emission measured with the Spitzer Space Telescope is not originating in a disk 
but likely in a jet/ouflow. 



- 32 - 



4 
3 

>^ 2 



X 



1 


-1 
0.8 
0.6 

^ 0.4 



X 



^ 1 

- 


i TWHya ; 

1 ^ 






■ , , , , 1 , , , , 


1 , , , , 1 , , , , - 





-50 

Av 



50 
heuo [km/s] 



-50 50 
[km/s] 



Fig. 4. — Comparison between observed and predicted (red solid line) line profiles for the 
group of transition disks (velocity in the stellocentric frame). Predicted line profiles are 
for the standard photoevaporative disk wind assuming an ionizing flux $ of 10^^ photons/s 
and the stellar mass and disk inclination in Table HI Models have been degraded to the 
spectral resolution of VISIR and scaled by the factors reported in the last column of Table H] 
to match the observed and predicted [Ne ll] line luminosities. The black thick line below 
the [Ne ll] emission shows the total uncertainty in the observed [Ne ll] peak position. Note 
the excellent match between the observed and predicted line profile from the almost face-on 
disk of TW Hya suggesting that the [Ne ll] emission originates from a photoevaporative disk 
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Table 1. Stellar properties and [Ne ll] fluxes measured with the Spitzer Space Telescope 



Source 


2MASS J^' 




Age 


Ref. 


[Nell] flux 


Ref. 






(km/s) 


(Myr) 




(xlO~i* erg cm~2 


s-i) 


HD 34700Aa* 


05194387+0539406 


68.0±0.2 


1-10 


1 




11 


TW Hya 


11015191-3442170 


12.2±0.5 


~10 


2,3 


4.5±0.3,5.7 


11,12 


CS Cha* 


11022491-7733357 


14.9±2 


~2 


4,5,6 


3.4±0.1,4.3 


11,13 


VW Cha* 


11080148-7742288 


17.2±0.5 


~2 


5,6 


2.8±0.4 


14 


T Cha 


11571348-7921313 


14.0±1.3 


>12 


5,7 


3.2±0.2 


14 


Sz 73 


15475693-3514346 


-3.3±2.5 


2.6-5.4 


8,9 


1.6±0.2 


14 


Sz 102 


16082972-3903110 


5 




10,9 


3.6±0.1 


14 


=>-The 2MASS 


source name includes 


the J2000 


sexagesimal, equatorial position in 


the form: hh- 



2003h 



mmssss-|-ddmmsss ijCutri et al 

'^Stellar heliocentric radial velocity 



* These stars are part of multiple ste llar systems. HD 34700 A is a spectroscopic binary and has two 
visual c ompanions iSterzik et al.l2005 ). CS Cha was recently found to be a long-period spectro scopic 

ahlll 



binary llGuenther et al. 



2007 ). VW Cha has a visual companion at ~0.7" i Brandner et 



19961) 
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Table 2. Summary of the observations. The observing time (in U.T.), the airmass, and 
the hehocentric radial velocity corrections (vheiio) are given at the beginning and at the end 

of the observations. 



Source 


Date 


U.T. 


texp 


Airmass 


Vheiio 


Calibrator 


texp 


Airmass 




yyyy-mm-dd 


hh:mm 


(s) 




(km/s) 




(s) 




TW Hya 


2008-02-20 


02:37/04:07 


3600 


1.3/1.1 


13.16/13.03 


HD 90957(K5III) 


240 


1.0(F) 


CS Cha 


2008-02-20 


05:00/06:29 


3600 


1.7/1.7 


12.01/11.97 


BV441(M0III) 






VW Cha 


2008-02-20 


07:07/07:15 


3480 


1.7/1.9 


12.06/12.03 


HD 92682(K3n) 


240 


1.6(P) 


VW Cha 


2008-02-20 


07:35/09:01 


3480 


1.7/1.9 


12.06/12.03 


HD 105.340(K2n) 


240 


1.7(F) 


HD 34700A 


2008-02-21 


01:22/02:51 


3600 


1.2/1.6 


-27.40/-27.55 


EPS Tau{G9.5III) 


240 


1.6(P) 


T Cha 


2008-02-21 


03:34/06:37 


7200 


1.9/1.7 


12.99/12.94 


BV441(M0in) 


240 


1.8(P) 














HD 105340(K2II) 


360 


1.6(F) 


Sz 102 


2008-02-21 


07:41/09:11 


3600 


1.3/1.1 


28.74/28.64 


HD 136422(K5ni) 


360 


1.2(P) 














HD 139127 


600 


1.1(F) 


CS Cha 


2008-02-22 


02:23/04:08 


4140 


1.9/1.7 


12.03/12.0 


BV441(M0ni) 


240 


1.9(P) 


CS Cha 


2008-02-22 


04:28/06:00 


3600 


1.7/1.7 


11.99/11.95 


HD 92682(K3n) 


360 


1.6(F) 


Sz 73 


2008-02-22 


06:58/09:17 


5400 


1.4/1.0 


29.36/29.19 


HD 136422(K5ni) 


360 


1.4(P) 














HD 139127(K4.5ni) 


600 


1.1(F) 
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Table 3. Summary of the VISIR results. For [Ne ll] detections we assume a Gaussian for 
the hne profile and a first order polynomial for the continuum. The last column gives the 

velocity of the peak emission in the stellocentric frame. 



Source 


[Ne II] 


Line Flux* 


FWHM 


^peak 




detection? 


( X 10"^"^ erg cm~^ s~^) 


(km/s) 


(km/s) 


HD 34700A 


N 


< 0.5 






TW Hya 


Y 


3.8±0.3 


14.6±0.7 


-6.2±0.3 


CS Cha 




2.3±0.2 


27±2 


-3.3±0.7 


VW Cha 


N 


< 0.6 






T Cha 


Y 


2.2±0.3 


42±4 


-4±2 


Sz 73 


Y 


1.1±0.2 


60±8 


-99±3 


Sz 102 


N 


< 0.5 







*In the case of non-detections we provide 3a upper limits to the line flux 
assuming a width of 10 km/s, equal to the instrument resolution. If the 
line is broader the flux upper limit would increase proportionally with the 
hne FWHM by the factor (FWHM/10). 



^In this case the continuum emission is not detected 
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Table 4. Stellar and disk parameters used to generate the [Ne ll] profiles in Fig. HI Stellar 
ionizing fiux $ were fixed to 10^^ photons/s. The last column provides the ratio between 
the observed and predicted [Ne ll] line luminosity, which is also the scaling factor for the 
ionizing fiux $. The comparison between models and observations yields $ = 2.5, 13, and 
1.9x10^^ photons/s for TW Hya, CS Cha, and T Cha respectively with an unceratinty of a 

factor of 2-3. 



Source Disk inclination Distance Ref. Lobs/Lmod 

(Mo) (°) (pc) 



TW Hya 


0.7 


4 


51 


1,2,3 


2.5 


CS Cha 


0.9 


45* 


160 


4,5 


13 


T Cha 


1.5 


75 


66 


6,7 


1.9 



* Chosen to provide a good match to the [Ne ll] line centroid and 
FWHM 
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